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ABSTRACT 

The model is presented of H 2 O maser in the lio — loi line within the first excited 
vibrational state of the molecule around oxygen-rich asymptotic giant branch stars. 
It is suggested that the maser cloud is located in the inner layers of the circumstellar 
envelope where intense dust formation takes place. The calculations took into account 
rotational levels belonging to the five lowest vibrational states of the H 2 O molecule. 
The model predicts the gain values of the 658-GHz maser about 10“^"^ — 10“^^ cm“^ 
at H 2 molecule concentrations 10® — 10^^ cm“^ and at high ortho-H 20 concentrations 
> 10® cm“®. The gas temperatures 1000 - 1500 K are considered to be a necessary 
condition for the effective maser operation. Results are presented for other maser 
transitions of the excited vibrational states of the molecule. 
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1 INTRODUCTION 

Intermediate- and low-mass stars (with main sequence mass 
1-6 M0) cease their life on the asymptotic giant branch 
(AGB), where they combine a high luminosity (i.e. 3000 L© 
and higher) an d a low effective temperature, below 3000 K 
llHabindll996l ). Stars with higher mass become red super¬ 
giants and have ten times the size of AGB stars. Late-type 
stars lose a substantial fraction of their initial mass by stel¬ 
lar winds and are surrounded by extended envelopes. Masers 
provide an unique probe of structure and kinematics of the 
circumstellar environment of late-type stars. Maser emission 
from the three most common maser molecules - SiO, H2O, 
and OH - traces regions of the circumstellar envelopes of 
oxygen-rich evolved stars on scales from a few to seve ral 
hundred astronomical units, hereafter AU (lAlcoleall2004l ). 

Intense H2O maser emission in the 22.2-GHz 616 —^ 623 
transitioik] belonging to the ground vibrational state of the 
molecule has been found towards many hundreds of cir¬ 
cumstellar envelopes of Mira and s emi-regular long-period 
variable stars and r ed supergiants (IValdettaro et al]|200ll : 
[Richards et al.l [20121 ). Radio interferometry methods allow 

* E-mail:alex-nlO@yandex.ru 

^ Energy levels of the asymmetric rotor molecule H 2 O are de¬ 
noted JKaKc where J is the total rotational quantum number, Aa 
and Ac represent its projections along the molecular axes. The 
vibrational states are denoted by quantum numbers (:^i, 1 ^ 2 , r' 3 ), 
where u\ denotes symmetric stretch, V 2 - bending, and 1/3 - an¬ 
tisymmetric stretch, respectively. 


us to image H2O masers in this line at high resolution and 
the structure of the circumstellar envelopes can be studied in 
detail. In addition to the 22.2-GHz line, maser emission from 
a number of other ortho- and para-H20 (sub)millimetre 
wavelength transitions belonging to the gr ound and first ex¬ 
cited vibrational st ates has been detected llHumphrevsl2007l : 
iHunter et al.ll2007l) . The maser emission in rotational lines 
within the excited vibrational states is of particular interest. 
Given the excitation conditions, these lines must arise from 
the innermost hot regions of circumstellar envelopes. 

The emission from two rotational lines within the first 
excited vibrational state (010) of H2O was discovered by 
iMenten &: Melnik ll 19891 ) using the Institut de Radioas- 
tronomie Millimetrique (IRAM) telescope. The para-H2 0 
transition 440 —>■ 633 at a frequency of 96 GHz and ortho- 
H2O transition 650 —f 643 at 233 GHz were detected towards 
the supergiant VY CMa, the latter transition was also de¬ 
tected towards the semiregular variable W Hya. The emis¬ 
sion at 96 GHz towards VY CMa showed signs of maser 
action. 

iMenten fc Youn^ lll995l) detected strong H2O maser 
emission at a frequency of 658 GHz towards a sample of 10 
oxygen-rich evolved stars using the 10.4-m telescope of the 
Caltech Submillimeter Observatory. The line corresponds to 
the lio — loi rotational transition within the (010) excited 
vibrational state of the ortho-H2 0 molecule. The 658-GHz 
spectra show contiguous, broad, emission with asymmetric 
line shapes centred at the stellar velocity. The photon lu¬ 
minosities of the 658-GHz masers, assuming the isotropic 
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emission, ranged from 5 x 10^® to 3 x 10“^® s“^. For the de¬ 
tected sources the photon luminosity of the 658-GHz line 
exceeded that of the 22.2-GHz line and the SiO maser lines. 
iHunter et al.l (120071 ') reported the observations of the 658- 
GHz circumstellar masers obtained by the Submillimeter 
Array. They presented new detections from six Mira vari¬ 
ables including IK Tau. [Richards et al.l (l2014h presented At¬ 
acama Large Millimeter/submillimeter Array (ALMA) ob¬ 
servations of the 658-GHz water maser towards VY CMa. 
The interferometric observations with high angular resolu¬ 
tion enabled us to locate water masers with milliarcsecond 
precision. The 658-GHz masers have a complex distribution 
close to the star and are located mostly outside the outer rim 
of SiO masers in this source. The 658-GHz maser emission 
is more concentrated towards the stellar position than the 
emission in th e other H 2 O maser lin es at 22.2 GHz, 321 GHz 
and 325 GHz dRichards et al.l[2014l L 


Other lines of the (010) excited vibrati o nal st ate of the 
H 2 O molecule were sought. iFeldman et al.l (Il993ll detected 
the 523 —>■ 6 i 6 rotational line of ortho-H 2 0 at 336 GHz 
toward s VY CMa u s ing t he James Clerk Maxwell Tele¬ 
scope. iMenten et al.l ll2006li reported the detection of the 
emission in the ortho-H 20 lines 651 — 752 at 294 GHz 
and 523 —>■ 616 at 336 GHz towards VY CMa using 
the Atacama Pathfinder Experiment (APEX) telescope. 
Searches have been performed for three transitions with neg- 


ative results: para-H 2 0 transition 650 
I Menten et al.l 


s^para^ 

_ ail 20^. o rtho-H 2 0 423 - 

i Mvers fc Barre tiT 1982h. and orth o-H?0 4 i 4 
67.8 GHz dPetuchowski fc Bennettlll99lh . 


753 at 297 GHz 
330 at 12.0 GHz 
321 at 


IPeguchil dl977l 'l included the vibrationally excited state 
( 010 ) in calculations of the pumping of H 2 O masers in cir¬ 
cumstellar envelopes of late-type stars. His model predicted 
maser emission in 12.0-GHz, 67.8-GHz and other lines be¬ 
longin g to this excited vibrational state. lAlcolea fc Ment^ 
dl993tl predicted maser action in several H 2 O transitions 
belonging to the first excited vibrational state whose up¬ 
per level quantum numbers are = J and = 0 
or 1. However, the maser action in t he 658-GHz lin e was 
not predicted in t he cal culations of IPeguchil dl977l ~l and 
lAlcolea fc Menteiil dl993ll . This maser line was also not 
mentioned in the study of c ircum stellar H 2 O excitation by 
iPeguchi fc Nguven-Q-RieuI dl990ll . 


Until recently the modelling of the vibrational ex¬ 
citation of H 2 O molecule was hampered by the lack 
of collisional excitation rates for transitions between the 
ground vibrational sta te and excited vibrational states. 
iFaure fc JosselinI d2008ll presented collisional rate coeffi¬ 
cients for the (de)excitation by H 2 of the lowest 824 ro- 
vibrational levels of the H 2 O molecule in the temperature 
range 200 — 5000 K. These data can be used for modelling 
the water molecule excitation in warm astrophysical media. 


In this paper we study the physical conditions neces¬ 
sary for the generation of strong maser emission in the vi¬ 
brationally excited H 2 O line at 658 GHz in the circumstellar 
envelopes around oxygen-rich AGB stars. In our model the 
658-GHz maser is located in the dust formation zone in the 
inner layers of the circumstellar envelope. 


2 PHYSICAL PARAMETERS OF THE MODEL 


2.1 Inner region of circumstellar envelopes 

The star’s physical parameters adopted in our calculations 
correspond to those of the AGB star IK Tau (Table [1}. IK 
Tau is an M-type sta r and is at a distance of about 250 pc 
llOlofsson et al.|[l998l l. The an gular diameter of the stellar 
photosphere was estimated by iMonnier et aP ll2004h to be 
20 mas based on the combined Keck aperture-masking data 
and Infrared Optical Telescope Array (IOTA) data. 

IK Tau exhibits circumstellar maser emission in SiO, 
H 2 O, and OH molecules. The SiO maser emission is found 
to arise from a clumpy ring within two to three stel¬ 
lar radii from the star centre dBoboltz fc Piamondl l2005l : 
iMatsumoto et aLll2008l L The emission of the vibrationally 
excited H 2 O masers and SiO masers arises from a sim¬ 
ilar range of velocities, which c an be interpreted as evi¬ 
dence for a close physi c al origin ([ Menten fc MelnickI Il99ll : 
iMenten fc Youn^ Il995l : iHunter et al. 2003). The observa¬ 
tions of the variability of the SiO maser distribution re¬ 
veal that the medium in the inner layers of circumstellar 
envelopes is strongly turbulent. The tur bulent veloc i ties of 
the g as can be as high as several km s’ (It iil l 19861; IPe^ 

1201:^ . 

Dynamical models of the AGB star winds pre¬ 
dict shocks and large velocity gradients up to several 
km s“^ AU~^ i n the inner layers of circumstellar en¬ 
velopes dWoitk e 20061: jNowotnv. Hofner fc Aringeil I 2 OI 0 I : 


[Ireland. Scholz fc WoodI 2011 1. The gas temperature in 
these layers tends to be close to the radiative equilibrium 
temperatu re, except for some narrow, hot post-sh ock cool¬ 
ing zones dSchirrmacher, Woitke fc Sedlmavill2003l l. 

Hydrogen is contained mainly in the form of H 2 
molecules in the photosp here of cool late-type stars 
dGlassgold fc HugginsI Il983l l . Hydrogen molecules are de¬ 
stroyed by the dissociative shock waves propagating through 
the stellar atmosphere and are then re-formed on dust 
grains. For the sake of simplicity, the hydrogen is considered 
to be molecular in the calculations. The relative abundance 
of H 2 O molecules in the stellar wind depends on the carbon- 
to-oxygen rati o in the stellar ph otosphere and can vary over 
a wide range dCherchnefill200B l . 


2.2 Dust model 

The mass loss observed in evolved AGB stars is usually at¬ 
tributed to a two-stage process: atmospheric levitation by 
pulsation-induced shock waves followed by radiative accel¬ 
eration of dust grains, which tran sfer momentum to the sur¬ 
rounding gas through collisions dBladh fc Hofneill2012l L In 
order for a dust species to trigger a wind it has to form 
close to the stellar surface, within reach of the shock waves. 
Radiation-hydrodynamical models of winds of oxygen-rich 
M-type AGB stars suggest that the winds can be driven 
by pho ton scattering on Fe-free silicate gra ins of sizes 0.1- 
1 qm (iHofnerl I2OO8I : Ifiladh fc Hofnerl r2()12h . The observa¬ 
tions of three M-type AGB stars using the multi-wavelength 
aperture-masking polarimetric interferometry pointed to the 
presence of dust shells with large grains (about 0.3 qm in ra - 
dius) at radii of about two stellar radii (iNorris et al.ll20i^ L 
The inner boundary of the dust formation zone appears to 
be located in or just outside the SiO maser zone. 
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Model of the vibrationally excited H 2 0 maser at 658 GHz 3 


Table 1. Physical parameters of the model 


Stellar radius 

R, = 2.5 AU 

Stellar temperature^ 

T* = 2300 K 

Distance from star centre^ 

D = 12.5 AU 

Cloud thickness^ 

H = 1.25 AU 

Number density of H 2 

3 X 10® sg sS 10 “ 

Number density of ortho-H 2 0 

lO-i ^ AhjO ^ 2 X 10® 

Gas temperature 

600 < Tg ^ 1600 K 

Microturbulent velocity 

r^turb = 3 km 

Velocity gradient 

fcv = 5 km AU“^ 

Dust-gas mass ratio 

0 /d s: 0.005 


^ The effective tem perature of the photosphere of IK Tau 
dMonnier et al.|[2004h . ^ See Fig. 1. 


In our model, the forsterite MR 2 Si 04 was considered as 
the dust grain material l|Bladh fc Hofneril2012h : some calcu¬ 
lations were done for magnesium-iron silic ates. Here we used 
the refractive index d ata for dust from djager et al. I I 2 OO 3 I: 
IPorschner et al.lllQQRh 0. In onr calculations, the dnst par¬ 
ticle radius a was assumed to be 0.3 pm. The cross sections 
for the absorption of radiation by dust particles were calcu¬ 
lated using the Mie scattering theory. We used the numeri¬ 
cal co de published in the monograph bylBohren fc HuffmarJ 
dl983ll and modified by iDrain i d2004(i R The dust grains 
grow until the elements contribnting to the dust material 
are consumed due to the relative element abundances in the 
stellar atmosphere. The first element in Mg 2 Si 04 that will 
be completely consumed is Mg for a solar abundance of ele¬ 
ments. The value of the dust-gas mass ratio /a determined 
by the relative abundance of Mg is approximately 0.002. 
The parameter /a was varied over the range 0 - 0.005 in the 
calculations. 

The dust heating mechanism in the inner layers of the 
circumstellar envelopes is the absorption of stellar radiation. 
The dust loses thermal energy through emission of radiation. 
The dust emissivity £c(i^) was calculated in accordance with 
Kirchhoff’s radiation law; the function £c(i^) strongly de¬ 
pends on dust temperature Ta. The dust temperature Ta at 
distance D from the stellar centre was calculated from the 
solution of the equation of thermal balance: 

00 00 

47r J dv ec{v) = J du (1) 

0 0 

where Kc(i^) is the dust absorption coefficient, Itiy) is the 
intensity of stellar radiation, is the solid angle that cuts 
out the stellar disk on the celestial sphere at distance D. 
Here, an optically thin stellar atmosphere is assumed, where 
the incident intensity on the grains is direct star light. 


^ http://www.astro.uni-jena.de/Laboratory/Database/databases 
® http://code.google.com/p/scatterlib/wiki/Spheres 


Gas-dust 

cloud 
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Figure 1. Gas-dust cloud model. The distance from the maser 
cloud to the star centre is equal to D = 5i?t. 

2.3 Maser cloud 

There is no interferometric data for the (sub)millimetre H 2 O 
masers with suc h resolution that ind ividual maser clouds 
can be resolved llRichards et al.lf2014ll . We suggest that vi¬ 
brationally excited H 2 O maser radiation originates in com¬ 
pact gas clouds in the warm inner layers of circumstellar 
envelopes. The maser cloud dimensions were set to be equal 
to 0.1 of the distance from the parent star. Such cloud sizes 
were revealed for the 22.2-G Hz H 2 O masers base d on radio 
interferometry observations llRichards et al.l[2012ll . 

We consider the one-dimensional model of a flat gas- 
dust cloud (see Fig.[T]). The z coordinate axis passes through 
the star centre and is perpendicular to the cloud plane. 
The cloud sizes along two coordinate axes are much larger 
than those along the third 2 coordinate axis. We assume 
that there is a gas velocity gradient along the z axis, fcv = 
dv/dz > 0. The cloud consists of a mixture of H 2 and H 2 O 
molecules. He atoms, and dust particles and is in the radia¬ 
tion field of the parent star. For simplification, the physical 
parameters of the cloud (the gas and dust temperatures, 
the number densities of atoms and molecules, and the dust 
abundance) were assumed to be independent of the coordi¬ 
nates. However, the H 2 O level populations were considered 
as functions of the z coordinate. We assume that H « D 
(plane geometry). 

The gas temperature in the maser clouds must sat¬ 
isfy Tg < 1600 K, because there is dissociation of 

H 2 O molecules at higher temper atures ll Jeong et al.l 1200^ : 
iNesterenok fc Varshalovichl I2OIII ). The model parameters 
adopted in our calculations are given in the Table [T] 


3 CALCULATION OF THE H 2 O LEVEL 
POPULATIONS 

3.1 System of master equations for the molecule 
level populations 

In the stationary case, the system of equations for the level 
populations is 
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M 

^ ^ H“ Chi') 

k = l, k^i 

M 

-ni{z) ^ {Rikiz) + Cik) = 0, i= 1, ...,M - 1, (2) 

k = l, k^i 

M 

= 1, 

i=l 

where rii is the normalized population of level i, M is the 
total number of levels, Rik{z) is the rate coefficient for the 
transition from level i to level k through radiative processes, 
and Cik is the rate coefficient for the transition from level i 
to level k through collisional processes. The rate coefficients 
for radiative transitions Rik{z) are 

~ ^ikJiki^z') T ^ 

( 3 ) 

^ik^A RikJiki^z') ^ Si ^ 

where Aik and Bik are the Einstein coefficients for sponta¬ 
neous and stimulated emission, respectively; Si and Sk are 
the energies of the levels i and k, respectively; Jik{z) is the 
radiation intensity averaged over the direction and over the 
line profile. For the plane parallel geometry: 

OO 1 

Jikiz) = ^ j j (f)ik{z,^,v)I[z,^,v). (4) 

— OO —1 

Here, is the intensity of radiation at a frequency 

V in direction p = cosO, where 9 is the angle between the z 
axis and the radiation direction, (j)ik[z, ^,v) is the normal¬ 
ized spectral line profile. The fine intensity I{z,fj.,u) is the 
solution of equation of radiative transfer: 

dl{z,u,v) , . . , , , 

M-= -n{z,ij.,v)l{z,ij,,v) +e{z,ij.,v), (5) 

where e{z, n, v) is the emission coefficient, and k[z, p, v) is 
the absorption coefficient. Each of the coefficients e{z,^i,v) 
and k{z, fi, v) is the sum of the emission or absorption coef¬ 
ficient of the dust and the emission or absorption coefficient 
in a spectral line, respectively: 

£{z,^i,u) = £c{v) + ^AikNni{z)(l>ik{z,fj,,iy), 

47r 

K{z,fl,lz) = lic{v)+ (g) 

+ ^AikN ( —nk{z) - ni(2) ) (j>ik{z, n, v). 

OTT \gk J 

Here, N is the molecule number density; gi and gk are the 
statistical weights of the levels; A is the radiation wave¬ 
length. In these formulas, it is implied that level i lies above 
level k in energy, i.e., Si > Sk- The spectral profile of the 
emission and absorption coefficients in the laboratory frame 
of reference is: 

(l)ik{z, g,,u) = ^ik{iz - g,Uikv{z)/c) (7) 

where <j)ik{iz) is the normalized spectral line profile in the 
co-moving frame of the gas, Vik is the transition frequency, 
v{z) is the gas velocity along the a axis, here v{z) = k^z, k^ 
is a constant coefficient equal to the gas velocity gradient in 
the cloud. 


The normalized spectral line profile in the co-moving 
frame of the gas is: 


4>ik{v) 



X= {v - Vik)lC^Vik 


(8) 


where is the line profile width. The line profile width is 
determined by the spread in thermal velocities of molecules 
and turbulent velocities in the gas-dust cloud: 


A 2 2 2 

^^ik — ^ik -5 = Ux + '^turb? 

C 


( 9 ) 


where vt = \j 2 kT^lm is the most probable thermal veloc¬ 
ity of the molecules, k is the Boltzmann constant, m is the 
mass of the molecule, and Uturb is the characteristic turbu¬ 
lent velocity in the cloud. 

Let fj,c to be the cosine of the angle between the 2 ; axis 
and the radiation direction from the edge of the stellar disk. 
Ate = a/I - Ri/D^. The boundary condition for eq. © at 
2 = 0 is 1(0, At, tz) = /»(tz). At ^ fJ-c, and /(0,Ai,tz) = 0, 
Ate > A* ^ 0- The star is assumed to have a black-body 
radiation spectrum with photospheric temperature equal to 
2300 K. At 2 = we take the boundary condition to be 
I{H, g,u) = 0, g, < 0. 

There is no fine overlap between the transitions of the 
H 2 O molecule. The emission from inverted transitions is dis¬ 
regarded in our calculations of the energy level populations 
of the molecule. 


3.2 Spectroscopic data and collisional rate 
coefficients 


In our calculations, we took into account 411 rotational lev¬ 
els of ortho-H 2 0 and 411 rotational levels of para-H 2 0 be¬ 
longing to the five lowest vibrational states of the ground 
electronic state of the molecule. The first rotational levels 
of the excited vibrational states have the following energies: 
1594.8 cm“^ for the (010) vibrational state, 3151.6 cm“^ 
for (020), 3657.1 cm”i for (100), and 3755.9 cm"i for (001). 
The spectroscopic data for H 2 O molecule were ta ken from 
the HITRAN 2012 database (iRothman et al.ll2013li . The en¬ 


ergies of the rotational-vibrational transitions for the H 2 O 
levels under consideration are less than 4400 cm~'^ (2.3 pm) 
and correspond to the infrared and the radio band. Note 
that the peak of the intensity of black-body radiation Biy) 
at a temperature of T* = 2300 K is near 2.2 pm. 

The rate coefficients for collisional transitions between 
H 2 O levels i n coll isions of H 2 O with H 2 were taken from 
iFaure et al.l ll2007l') and iFaure fc Josselii] ll2008ll . The rate 
coefficients for collisional transitions between H 2 O lev¬ 
els in inel astic collisions of H 2 O with He a toms were 
take n from [Gree n. Maluendes fc McLear] (Il993h . The data 
from Icireen. Maluendes fc McLe^ ( 1993h contain the col¬ 
lisional rate coefficients for the lowest 45 levels of ortho- 
H 2 O molecule and 45 levels of para-H 2 0. The collisional 
rate coefficients for the transitions including higher H 2 O 


iGreen, Maluendes & McLeai 

J (I 1993 

proposed bv Faure & Josselir 

J ll2008h 


using the algorithm 
for the rate coefficients 
for H 2 O and H 2 collisions. In our calculations, we used the 
relax ation rates of excited vibrational states of H 2 O by He 
from iKung fc Centeil lll975h . The ratio of the summed colli¬ 
sional rate coefficients out of any H 2 O level for H 2 O - He and 
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H 2 O - H 2 collisions lies in the range 0.1 - 1. The abundance 
of helium atoms relative to the total number density of hy¬ 
drogen nuclei was assumed to be A^He/ 2 A^H 2 = 0.1. Thus, 
the contribution of H 2 O - He collisions to the collisional 
excitation of H 2 O molecule is about 10 per cent. 


3.3 Numerical calculations 

The system of master equations for the H 2 O level pop¬ 
ulations eq. © and radiative transfer equation in the 
medium, eq. form a system of non-linear equations 
for the level populations of the molecule. This system of 
equations was solved by th e accelerated A-iteration method 
dRvbicki fc Hummeilll99ll i. The Feautrier method was used 
for solving th e radiative transfer equation (iFeautrierj 1 19641 : 
|Peraiahll2004l '). The collisional and radiative transitions be¬ 
tween ortho and para spin isomers of H 2 O molecule are for¬ 
bidden, and the level populations of the molecules are cal¬ 
culated independently. 

The cloud in our numerical model was broken down into 
layers parallel to the cloud plane. The molecular level pop¬ 
ulations within each layer were constant. The thickness of 
the near-surface layer was chosen to be A^s = cm. 

The optical depth for any H 2 O line and for any considered 
direction in the near-surface layer was less than 1. The thick¬ 
ness of each succeeding layer into the cloud was larger than 
that of the preceding one by a constant factor. The number 
of layers into which the cloud was broken down is 100. The 
calculations showed that the decrease of the surface layer 
thickness (by a factor of 10 ) or the increase of the number 
of layers lead to a small effect on the population inversion 
of the transitions in question (< 1 per cent). 

The intensity is averaged over opposite directions in the 
Feautrier calculational scheme; the range of values for the 
parameter p is [0; 1]. The number of angles used in angu¬ 
lar integration was set to be 15. The minimal value of the 
parameter /r for which the radiative transfer equation was 
solved is equal to 0.05, the corresponding dji is equal to 0.1. 
The maximal value of the parameter fj, is equal to 1 and cor¬ 
responding djj. = 1 — fic- Other angle values are distributed 
uniformly over the interval. The range of values for the pa¬ 
rameter X for each line (see eq. o was chosen to be [—5; 10]; 
the centre of the interval is shifted to positive values in or¬ 
der to take into account the large velocity gradients in the 
cloud (fcv ^ 0). The discretization step of the parameter x is 
0.2. The model was tested with the number of frequency or 
angle points twice as much as listed above; the changes in 
the population inversions were negligible (~ 0.1 per cent). 

An additional acceleration of the iterative series was 
achieved by applying t he convergence optimization method 
proposed bv iNd (Il974l ~). The convergence criterion for the 
iterative series was the condition on the maximum relative 
increment in level populations for two successive iterations, 

max|Ani/ni| < 10~"'. (10) 

i 

The error of the population inversion corresponding to this 
convergence criterion is < 1 per cent for the parameter val¬ 
ues of physical interest. 

In our calculations, we used the al gorithms for solving 
system s of linear equations published in R ybjcki fc Hummer! 
lll99lfl and in the book by IPress et al. I (1199711 . The calcu¬ 



Figure 2. The gain of the 658-GHz maser line as a function of 
the cloud depth. The results are given for three cases: a. The 
411 rotational levels belonging to the ground and four excited 
vibrational states of the molecule are taken into account (solid 
line), b. The 271 rotational levels belonging to the ground and 
the first excited vibrational states of the molecule are taken into 
account (dashed line), c. The 411 rotational levels belonging to 
the ground and four excited vibrational states of the molecule 
are taken into account, and the stellar radiation is not taken into 
account in the radiative transfer (dashed-dotted line). 


lations were performed on the supercomputer of the Saint 
Petersburg branch of the Joint Supercomputer Center of the 
Russian Academy of Science^ 


4 RESULTS 


4.1 Gain of the 658-GHz maser line 


The expression for the gain of the transition i ^ jf at the 
line centre in a direction along the cloud plane is 




X^AikN 




( 11 ) 


The cloud-averaged gain is calculated from the formula 



( 12 ) 


Inverted transitions with a small gain 7 < cm“^ 

were ignored in the result analysis. 

The dependence of the gain of the 658-GHz maser line 
on depth into the cloud is shown in the Fig.[2l The adopted 
parameters are Tg = 1100 K, Nh 2 = 5 x 10® cm“®, ortho- 
H 2 O concentration AH 2 O = 3 x 10® cm“® and dust-gas 
mass ratio fd = 0 . 002 , and the Mg 2 Si 04 dust material 


http://scc.ioffe.ru/ 
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Figure 3. The gain of the 658-GHz maser line as a function of 
the number density of hydrogen molecules and the number 
density of ortho-H 2 0 molecules A1 h 2 0- The contours of equal gain 
are shown; the gain in cm“^ is indicated near each curve. The 
inclined dash-dot lines correspond to the points for which the 
relative ortho-H 20 abundance xhjO = .^^H 2 o/ 2 AfH 2 equals 10“® 
and 10“^. The adopted parameters are Tg = 1100 K and = 
0 . 002 . 


was considered. The results are also presented for calcu¬ 
lations with smaller number of ro-vibrational levels of the 
molecule and for calculations without taking into account 
the stellar radiation in the radiative transfer equation. The 
differences between calculated results imply that the inclu¬ 
sion of the high-lying excited vibrational states is impor¬ 
tant in modelling of the emission from the lines belonging 
to the first ex cited vibrational state. Analogous conclusion s 
were made bvIConzalez- Alfonso. Neufel d fc MelnickI ll2007ll . 
iMaercker et al.l 1 2009f ). iDecin et al.l ( 201Clf l. who modelled 
the H 2 O line emission from the ground vibrational state 
of the molecule in the circumstellar envelopes of AGB 
stars. The vibrational excitation of the molecule takes place 
through collisions and through absorption of stellar radia¬ 
tion. The gain is substantially lower when the stellar radi¬ 
ation is taken into account; the cloud averaged gain 7 at 
the parameters in question is about 0.8 x 10 “^® cm“ 


diation, respectively (with a full set of ro-vibrational levels). 
Analogou s results were found for the 22.2-GHz circumstellar 
maser bv iNesterenokI (l2013lJ) . 

Fig. |3] demonstrates the calculated cloud-averaged gain 
7 of the 658-GHz maser line as a function of the number 
densities of hydrogen molecules H 2 and ortho-H 2 0 molecules 
in the cloud. The gain is 7 « 10“^® cm“^ at H 2 number 
densities about 3 x 10® to 5 x 10^® cm“® and ortho-H 20 
number densities about 3 x 10® to 10® cm“®. The excitation 
temperature of the maser levels is Texc ~ —10® K for this 
range of physical parameters. The length of the amplification 
region along the line of sight probably does not exceed, or 
is comparable to, the distance from the cloud to the star D. 


Figure 4. The gain of the 658-GHz maser line as a function of the 
dust-gas mass ratio. The adopted parameters are Tg = 1100 K, 
Ajjj = 5 x 10® cm“^ and ortho-H20 concentration A^HjO = 
3 x 10® cm“®. 


Thus, the optical depth in the maser line can be as high as 
r « 15 in the model in question. 

The maser becomes saturated if the transition rates 
between the maser levels are overwhelmingly determined 
by the induced radiative transitions in the maser line. As 
the degree of maser saturation increases, the maser lumi¬ 
nosity per unit gas volume a pproaches its limi ting value 
and the line profile broadens (IStrel’nitskiH llQTSh . The ob¬ 
served contiguous, broad emission of the 658-GHz masers 
may be interpreted as i mplying that the mase rs have a 
high degree of saturation (iMenten fc Younall995fl . The lim¬ 
iting m aser luminosity can be fou nd using the expression 
given in iNeufeld fc MelnickI lll99lll . In our model, the lim¬ 
iting photon luminosities of the maser he in the range 
10—100 cm“® s“® at the gas densities and water abundances, 
at which the gain values are high 7 > 5 x 10“^'^ cm“®. The 
limiting photon luminosity of the maser cloud is estimated 
Tiim ~ ~ 10 "*® s“® in the model in question (the 

strong maser emission is suggested to emanate from a cloud 
predominantly in one direction). We can estimate the beam¬ 
ing angle of a cylindrical maser AH ~ (H/D)^ ~ 10“® sr. 
The isotropic photon luminosity is Lis « 47 rLiim/AH ~ 
10"^“* s~® in the model in question. The calculated values of 
the photon luminosity can roughly reproduce the observed 
intensities of the maser radiation. 

Fig. 2] shows the dependence of the gain of the maser 
line 7 on dust-gas mass ratio. The results are presented 
for three dust materials. The radiative equilibrium tempera¬ 
tures for Mg 2 Si 04 , Mgo.yFeo.sSiOa and Mgo.sFeo.sSiOs dust 
grains are 340 K, 440 K and 660 K, respectively. The gain 
of the maser line increases with increasing dust abundance 
in the case of Mg 2 Si 04 and Mgo.rFeo.sSiOa dust material. 
For the Mgo.sFeo.sSiOs dust material the opposite is true 
- the gain decreases with increasing dust-gas mass ratio. 
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The radiation of warm dust suppresses the population in- 
ve rsion in the maser line. An alogous results were reported 
bv lYates. Field fc Gravl (119971 1 for collisionally pumped H 2 O 
masers from the ground vibrational state (22.2-, 325-GHz 
and other lines). The masers that are formed by collisional 
mechanism are quenched by the strong infrared continuum. 
This fact explains also the substantial decrease of the gain 
of the 658-GHz maser when the stellar radiation is taken 
into account in radiative transfer (see Fig. O- 

4.2 Gain of the other vibrationally excited masers 

The population inversion for other lines of H 2 O molecule 
belonging to the (010) vibrational excited state were con¬ 
sidered. The gain of the ortho-H 2 0 transition 4 i 4 —^ 821 at 
a frequency of 67.8 GHz was found to be 7 > 10“ cm“^ 
at ortho-H 20 concentration A^H 20 > 3 x 10® cm“® and 
Tg = 110 0 K. The maser emission in this line was predicted 
earlier bv lDeeuchil lll977ll . The gain of the ortho-H 20 tran¬ 
sitions 650 —>■ 643 at a frequency of 233 GHz and 651 — 752 
at a frequency of 294 GHz was found to be 7 > 10“^^ cm“^ 
at YH 2 O > 5 X 10® cm“®. There is no population inversion 
for the ortho-H 20 transition 623 —^ 616 at a fre quency of 
336 GHz in our model. Note. lMenten et al.l (l2006l l observed 
the 336-GHz line towards VY GMa and concluded, based on 
the simple considerations, that the line has a thermal excita¬ 
tion. The calculations showed that the population inversion 
for the transition lio — loi at a frequency of 794 GHz within 
the (020) vibrational state of the molecule can be high. The 
gain of the transition is found to be 7 > 10 “^^ cm“^ at the 
parameters Y’h 2 > 5 x 10® cm“®, Y’h 20 > 5 x 10® cm“® and 
Tg = 1100 K. This line is a twin of the 658-GHz maser line 
but is located within the higher excited vibrational state. 
The calculations did not predict any strong vibrationally ex¬ 
cited (sub)millimetre maser for the para-H 20 molecule. The 
gain of the para-H 20 vibrationally excited line 440 633 at 

a frequency of 96 GHz was found to be low 7 < 10“^'* cm“^. 

Fig. Opresents the dependence of the gain of the vibra¬ 
tionally excited maser lines on gas temperature. The gain 
depends strongly on gas temperature for the 67.8-GHz, 658- 
GHz and 794-GHz masers, while for other transitions the 
gain is almost constant. 

According to our calculations, there is weak depen¬ 
dence of the gain of the maser lines on the gas velocity 
gradient. The variation of the velocity gradient from 0 to 
9 km s“^ AU“^ leads to about 5 — 10 per cent change in the 
gain of the maser lines in question. The sensitivity of the 
calculated gain values to a change in the cloud height was 
investigated: the change in the cloud height by 10 per cent 
leads to a 1-5 per cent change in the calculated gain values. 


5 DISCUSSION 

The majority of the H 2 O masers are collisionally pumped 
dYates. Field fc GravlHOl^ . The collisional pumping cycle 
consists of collisional excitation of high-lying ro-vibrational 
levels of the molecule followed by radiative deexcitation 
of the levels (the emission of ’sink’ photons). The ’sink’ 
photons are absorbed by dust or escape from the reso¬ 
nance region to complete the pumping cycle. The collisional 
maser pumping is most efficient under conditions when the 



Figure 5. The gain of the maser lines as a function of the gas tem¬ 
perature. The line frequency in GHz is indicated near each curve. 
The adopted parameters are = 0.002, Yhj = 5 x 10® cm“®, 
ortho-H 2 0 concentration Y 0 H 2 O = 3 X 10® cm“® and para-H 2 0 
concentration YpH 20 = 10® cm“®. 


dust temperature is much lower than the gas temperature 
l[Bolgovaj_Stoelnitskii fc S hmeldlll977l : IChandra et al.l [19841 : 
lYates. Field fc Gra\^Il997l l. According to our calculations, 
the gain of the 658-GHz maser line rises signihcantly with 
increasing abundance of cold dust and with increasing gas 
temperature (see Figs |4] and [5]). These facts indicate that 
the 658-GHz maser line is collisionally pumped. There is not 
strong dependence of the gain on the gas temperature for the 
96-GHz, 233-GHz and 294-GHz maser lines. Moreover, the 
gain of these maser lines is substantially lower when the stel¬ 
lar radiation is not taken into account in radiative transfer 
equation. We suggest that these masers have an important 
radiative pumping component. 

According to the present calculations, gas densities 
about 10 ® — 10 ^® cm“® and relative H 2 O abundances 10 “® — 
10“'* are necessary for the effective pumping of the 658- 
GHz maser. Such values of the water abundance were con- 
sidered in the models of 22.2-GHz circ umstellar masers 
liHumphrevs et ahll200ll : lNesterenokll2013^ . The gas density 
in the stellar wind is about Nh 2 ~ 10 ® cm“® at distances 
D ~ 5 stellar radii from the star in a model of the circum¬ 
stellar envelope of oxygen-rich Mira variable with mass loss 
rate M rs 7 x 10“® M© year“* (Ijeong et al.ll200jl . The gas 
density and water abundance in the maser clouds can differ 
substantially from the values calcula ted in the framework 
of the isotropic c ircumstellar wind dRichards et al.l I2OI2I : 
lNesterenokll2()13^ ~) . 

The gas temperature in the maser cloud must be suf- 
hciently high for the effective collisional excitation of ro- 
vibrational levels of the H 2 O molecule. At the densities 
10 ® — 10 ** cm“®, the efficiency of the energy exchange be¬ 
tween gas and radiation field is at maximum which results 
in a close coupling of the gas to the condition of radiative 
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equilibrium dSchirrmacher, Woitke fc Sedlmavill2003h . The 
658-GHz maser emission in the circumstellar envelope of 
supergiant VY CMa is found to arise from an extended 
region at distances up to 40 stellar radii from the star 
llRichards et al.ll2014h . The radiative equilibrium tempera¬ 
tures are too low for the pumping of 658-GHz masers at 
such distances. Thus, the local heating mechanisms must 
play a significant role in the gas energetic equilibrium in the 
maser clou ds. These gas heating mechanisms can be H 2 - 
formation (ISchirrmacher. Woitke fc Sedlma^ 2003ll. dissi¬ 
patio n of turbulent motions of the gas ( Strelnitski et all 
|2002| ~). dissi pation of magnet o -hydr odynamic and acous¬ 
tic waves dPijpers fc HearrJ Il989li . The gas temper¬ 
ature in the post-shock layers, where the hydrogen 
molecule formation takes place, can be higher by 300 
- 400 K than the radiative equilibr ium temperature 
(ISchirrmacher. Woitke fc Sedlmavil 120031 '). Note that some 
maser emission in VY CMa originat es from elongated fea - 
tures characteristic of post-shock gas dRichards et al.ir2014h . 

Only those dust particles whose radiative equilibrium 
temperature is lower than the condensation temperature 

i which is about 1000 - 1100 K for silicates) can survive 
Gailll201^ 1. The strong candidates for dust material in the 
inner layers of circumstellar envelopes of oxygen-rich AGB 
stars are silicates with poor Fe abundance dBladh fc Hofnen 
I 2 OI 2 I) . Fe-free silicate material is more efficient at emitting 
the radiation than absorbing it. The relatively cold Fe-free 
silicate dust makes the collisional pumping of the masers 
efficient. 

The time variation of H 2 O circumstellar masers in 
22.2-GHz line is generally attributed to a change of 
the gas den sity and the velocity field of a circumstel¬ 
lar envelope ( Deg^hi| l 19771 : iRudnitskii fc Pashchenk3l2005l : 
IShintani et al.l 120081 '). We suggest that the variation of the 
infrared stellar luminosity can affect significantly the maser 
luminosities because of a strong coupling of the gas to the 
condition of radiative equilibrium. The lower the stellar lu¬ 
minosity, the lower the gas temperature and the less effective 
the collisional pumping of the masers. Moreover, the strong 
infrared radiation field affects the pumping mechanism and 
lowers the gain of the collisionally pumped masers (for ra- 
diatively pumped masers the effect has an opposite sign). It 
is believed that the inner radius of the H 2 O maser emission 
zone in a circumstellar envelope is determined by the dis¬ 
tance from the star where the gas density in the maser cloud 
falls below the collisional quenching density for the maser 
d Cooke fc Elitzurl 1 19851 : [Richards et al.l l2012l f. We suggest 
that the stellar radiation is another factor, that can affect 
significantly the maser pumping mechanism and quench the 
collisionally pumped masers in the inner layers of circum¬ 
stellar envelopes of late-type stars. 


6 CONCLUSIONS 

We investigate the physical conditions necessary for the gen¬ 
eration of the vibrationally excited H 2 O maser at 658 GHz in 
the inner layers of the circumstellar envelopes of AGB stars. 
It is shown that the high-lying excited vibrational states and 
the stellar radiation have to be included in the calculations 
of the level populations belonging to the first excited vibra¬ 
tional state. There is strong dependence of the gain of the 


658-GHz maser line on gas temperature and on dust compo¬ 
sition and abundance. The maser is found to be collisionally 
pumped. 

The gain of other maser transitions belonging to excited 
vibrational states were calculated. The ortho-H 2 0 masers 
at 67.8 GHz and 794 GHz were found to be collisionally 
pumped. The ortho-H20 masers at 233 GHz and 294 GHz 
and the para-H 2 0 maser at 96 GHz were found to have an 
important radiative pumping component. We suggest that 
the maser transition at 67.8 GHz might prove to be de¬ 
tectable from the ground, despite the fact that this maser 
was searched for earlier without success. 

It will be of particular interest to monitor observations 
of the stellar 658-GHz maser and other vibrationally excited 
maser lines and to determine of the time variability proper¬ 
ties of the H 2 O masers. 
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